Mon. Not. Ft. Astron. Soc. 000, [|-?? (2000) Printed 5 February 2008 (MN WT&L style file vl.4) 



A theoretical study of the mass temperature relation for 
clusters of galaxies 



(N ' 

O : A. Del Popolo, 1 ' 2 ' 3 
O ■ 
(N ■ 



> 
ON 

o 

O 

Oh 
6 



1 Dipartimento di Matematica, Universitd Statale di Bergamo, Piazza Rosate, 2-1 24129 Bergamo, ITALY 

2 Feza Giirsey Institute, P. 0. Box 6 Qengelkoy, Istanbul, Turkey 



a ■ 

Si 

\o ■ 

3 Bogazigi University, Physics Department, 80815 Bebek, Istanbul, Turkey 



Accepted ??? Received 2000 July 24; in original form 



ABSTRACT 

I derive the mass-temperature relation and its time evolution for clusters of galax- 
ies in different cosmologies by means of two different models. The first one is a modi- 
fication and improvement of a model by Del Popolo & Gambera (1999), namely based 
upon a modification of the top-hat model in order to take account of angular momen- 
tum acquisition by protostructures and of an external pressure term in the virial the- 
orem. The second one is based on the merging-halo formalism of Lacey & Cole (1993), 
accounting for the fact that massive clusters accrete matter quasi-continuously, and is 
an improvement of a model proposed by Voit (2000) (herafter V2000), again to take 
account of angular momentum acquisition by protostructures. The final result is that, 
in both models, the M-T relation shows a break at T ~ 3 — 4keV. The behavior of the 
M-T relation is as usual, M oc T 3 / 2 , at the high mass end, and M oc T 7 , with a value 
of 7 > 3/2 depending on the chosen cosmology. Larger values of 7 are related to open 
cosmologies, while ACDM cosmologies give results of the slope intermediate between 
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the flat case and the open case. The evolution of the M-T relation, for a given M v i r , 
is more modest both in flat and open universes in comparison to previous estimate 
found in literature, even more modest than what found by V2000. Moreover the time 
evolution is more rapid in models with L = than in models in which the angular 
momentum acquisition by protostructures is taken into account (L =/= 0). The effect 
of a non-zero cosmological constant is that of slightly increasing the evolution of the 
M-T relation with respect to open models with L =/= 0. The evolution is more rapid 
for larger values (in absolute value) of the spectral index, n. 

The mass-temperature relation, obtained using the quoted models, is also compared 
with the data by Finoguenov, Reiprich & Bohringer (2001) (hereafter FRB). The 
comparison shows that the FRB data is able to rule out very low fio models (< 0.3), 
particularly in the open case, and that better fit are obtained by ACDM models and 
by CDM models with Q > 0.3. 

Key words: cosmology: theory - large scale structure of Universe - galaxies: formation 

1 INTRODUCTION 

In the past decade, observations of clusters of galaxies (e.g, ROSAT, ASCA) have shown the existence of a tight correlation 
between the total gravitating mass of clusters, M to t, their X-ray luminosity (I/x), and temperature (Tx) of the intra-cluster 
medium (ICM) (David et al. 1993; Markevitch 1998; Horner, Mushotzky & Sharf 1999 (hereafter HMS)). The importance of 
these relations is due to the fact that cluster masses are difficult to measure directly, and when comparing cluster observations 
with models of structure formation a surrogate for cluster mass is usually used. Since M to t compares with the ICM temperature 
measurements which can be obtained through X-ray spectroscopy, this explains the importance of an M-T relation. Numerical 
simulations and theoretical studies have motivated the existence of such kind of relations, whose comparison can provide 
strong constraints on the prevailing cosmological models. 

For what concerns N-body simulations, in particular, Evrard, Metzler & Navarro (1996), and Eke, Navarro & Frenk 
(1997) showed that there is a tight relation between the mass of a cluster and its global X-ray temperature in gas-dynamical 
simulations, irrespective of the assumed cosmological model and of the state of the cluster. Other N-body simulations by 
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A theoretical study of the mass temperature relation for clusters of galaxies 3 

Mathiesen & Evrard (2001) have shown that the ICM X-ray temperatures show a tight (< 20% scatter) correlation with the 
total cluster mass components, correlation also seen in observations (Mohr, Mathiesen & Evrard 1999). 

From the theoretical point of view, within the framework of pure gravitational infall (e.g., Lilje 1992) by means of a 
top-hat model and simple arguments, based on the virialization density, suggest a self-similar relation, Mtot oc Tx Q confirmed 
by the results of N-body simulations (Evrard, Metzler & Navarro 1996; Eke et al. 1998). The proportionality coefficient of the 
law has been usually fixed by using numerical simulations to calibrate it (Evrard et al. 1996; Thomas et al. 2001). Afshordi 
& Cen (2001) (hereafter AC), through of energy conservation and assuming nearly spherical collapse, were able to obtain a 
normalization factor in agreement with hydro-simulations, but about 50% lower than X-ray mass estimates (see also V2000) 
(this discrepancy with observational normalizations is probably due to some effect, happening at small scales (AC), which 
the model does not take into account). For sake of completeness, it is important to stress that the results from different 
observational methods of mass measurement are not consistent with one another and with simulation results (e.g., HMS; 
Nevalainen, Markevitch & Forman 2000 (hereafter NMF), FRB). While mass estimates from galaxy velocity dispersion are 
consistent with simulations (HMS), X-ray mass estimates can differ by 80% from the prediction of hydro-simulations. 

Another drawback of the standard top-hat model has been stressed by Voit & Donahue (1998) (hereafter V98) and V2000. 
Using the merging-halo formalism of Lacey & Cole (1993), which accounts for the fact that massive clusters accrete matter 
quasi-continuously, they showed that the M-T relation evolves, with time, more modestly than what expected in previous 
models and this evolution is even more modest in open universes. 

Finally, there is another problem that, added to the previously quoted ones, urges for a better theoretical understanding 
of the M-T relation: recent studies have shown that the self-similarity in the M-T relation seems to break at some keV (NMF; 
Xu, Jin & Wu 2001). By means of ASCA data, using a small sample of 9 clusters (6 at 4 keV and 3 at ~ 1 keV), NMF has 
shown that Mtot oc j£- 79±0 - 14 f or the whole sample, and M to t oc T^/ 2 excluding the low-temperature clusters. Xu, Jin & Wu 
(2001) has found M to t oc T^ 60 ± 04 ( us i ng the /3 model), and M tot oc T^ 81 ± 014 by means of the Navarro, Frenk & White 
(1995) profile. FRB have investigated the T-M relation in the low-mass end finding that M oc T~ 2 , and M oc r~ 3 / 2 a t the 
high mass end. This behavior has been attributed to the effect of the formation redshift (FRB) (but see Mathiesen 2001 for a 
different point of view), or to cooling processes (Muanwong et al. 2001) and heating (Bialek, Evrard & Mohr 2000). AC has 
shown that non-sphericity introduces an asymmetric, mass dependent, scatter for the M-T relation altering its slope at the 
low mass end (T ~ 3 keV). 

* a = 3/2 is almost unrelated to the settings of the cosmological parameters 
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In this paper, I derive the mass-temperature relation, and its time evolution, for clusters of galaxies, in different cosmologies by 
means of two different models: the first one is a modification and improvement of a model by Del Popolo & Gambera (1999), 
based upon a modification of the top-hat model in order to take account of angular momentum acquisition by protostructures 
and of an external pressure term in the virial theorem. The second one is an improvement of a model proposed by V2000, 
again to take account of angular momentum acquisition by protostructures. Both models show that the M-T relation is not 
self-similar. A break is present in the quoted relation at T ~ 3keV and at the lower mass end the power law index of the M-T 
relation is larger than a = 3/2 even in flat universes. The slope of the power-law index depends on the considered cosmology. 
The two models also agree in predicting a more modest time evolution of the quoted relation in comparison with the results 
of previous models, and again depending on the cosmology. 

The plan of the paper is the following: in section (2) I describe the top-hat modified model (hereafter THM) and 
in section (3) the modified V2000 model (hereafter VM). Section (4) describes the results and section (5) is devoted to 
conclusions. 

2 REVIEW OF THE TOP-HAT MODEL 

As previously quoted, numerical methods and simple scaling arguments suggest that the X-ray temperature of clusters, Tx, 
can be directly related to their masses. In fact, assuming that the mean density within a radius r is A times the critical 
density, pb, one can express the mass in the given radius as: M(A) oc p^ 1 A -1 / 2 . Defining the quantity A v i r as the 
density contrast of a spherical top-hat perturbation just after collapse and virialization, in the "vein" of the Press-Schechter 

q 1 2 1/2 1/2 

approach, the mass can be expressed as M v j r = M(A v i r ) and so M v i r oc T x p h A vir , which gives the quoted relation 
between mass and temperature. 

In the following, I'll use Del Popolo & Gambera (1999) in order to get the M-T relation. The model is fundamentally a 
modification of the top-hat model in order to take account of angular momentum acquisition by protostructures and uses a 
modified version of the virial theorem in order to include a surface pressure term (V2000, AC). This correction is due to the 
fact that at the virial radius r v i r the density is non-zero and this requires a surface pressure term to be included in the virial 
theorem (Carlberg, Yee & Ellingson 1997) (the existence of this confining pressure is usually not accounted for in the top-hat 
collapse model). 

The equation governing the collapse of a density perturbation taking account angular momentum acquisition by proto- 
structures can be obtained using a model due to Peebles (Peebles 1993) (see also Del Popolo & Gambera 1998, 1999). 
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A theoretical study of the mass temperature relation for clusters of galaxies 5 

Let's consider an ensemble of gravitationally growing mass concentrations and suppose that the material in each system 
collects within the same potential well with inward pointing acceleration given by g(r) (see Del Popolo & Gambera 1998). 
We indicate with dP = f(L,rv r ,t)dLdv r dr the probability that a particle can be found in the proper radius range r, r + dr, 
in the radial velocity range v r — f, v r + dv r and with angular momentum L = rvg in the range dL. The radial acceleration 
of the particle is: 

dv r _ L 2 (r) _ L 2 (r) GM 

dt M 2 r 3 9(T) M 2 r 3 r 2 ( ' 

Eq. (jl|) can be derived from a potential and then from Liouville's theorem it follows that the distribution function, /, satisfies 
the collisionless Boltzmann equation: 



dt dr dv r 



Li . . 



(2) 



Assuming a non-zero cosmological constant Eq. (|l|) becomes: 
dv r GM L 2 (r) A 

~dt ~ ~ + MV + I T ( ' 

(Peebles 1993; Bartlett & Silk 1993; Lahav 1991; Del Popolo & Gambera 1998, 1999). Integrating Eq. (§) we have: 
1 fdr\ 2 GM f L 2 A 2 

su) = —+jM^ dr+ r +e (4) 

where the value of the specific binding energy of the shell, e, can be obtained using the condition for turn-around, 4? = 0. 

The CDM spectrum used in this paper is that of Bardeen et al. (1986) normalized to reproduce the observed abundance of 
rich cluster of galaxies (e.g., Bahcal & Fan 1998) Filtering the quoted spectrum on clusters scales, Rf = 2>h~ x Mpc, I obtained 
the total angular momentum acquired during expansion, as described in Del Popolo & Gambera (1998, 1999) (more hints on 
the model and some of the model limits can be found in Del Popolo, Ercan & Gambera (2001)). 

The temperature-mass relation can be obtained using the virial theorem, energy conservation and using Eq. (^). In Eq. 
(^), three forms of potential are present and the virial theorem can be written in the form: 

(K) = -±(U G ) - (U h ) + (U A ) (5) 

(Landau & Lifshitz 1960; Lahav et al. 1991, Del Popolo & gambera 1999), where, using Bartlett & Silk (1993) notation, the 
meaning of () is the time-average of the quantity considered, Ug the gravitational potential energy, Ul that connected to 
the angular momentum and Ua that due to A. For example, in the case of the gravitational potential energy, we can write 
Ug — — r^> wnere r eff is the time-averaged radius of a mass shell. As previously mentioned, when a surface pressure term 
is added, Eq. (|^) must be modified. Assuming that: 
© 2000 RAS, MNRAS 000, 
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Fig. (1) Fig. (2) 

Figure 1. M-T relation predicted by the modified top- hat model. The solid line is the prediction of Eq. (^), for f2yv = 0, f2rj = li shifted 
downwards, similarly to AC, to fit the FRB observational data in the massive end. The dotted and dashed line represents the prediction 
of the same equation for L ^ 0, VIa = 0, VLq = 1 and L ^ 0, Qa = and Qq = 0.3, respectively. 

Figure 2. As Fig. 1 but now Qjy 0. The solid line is equivalent to the solid one in Fig. 1. The dashed line represents the case of a flat 
universe with L ^ 0, A ^ 0, (Q = 0.3, Q A = 0.7). 



(K) + (E) = 3P ext F 



(6) 



(AC), where P ex t is the pressure of the outer boundary of virialized region and V is the volume, and by using the same 
assumption of AC, namely: 



3P oxt V = -uU 



(7) 



where v is a coefficient and U is the total potential. With the introduction of the surface term Eq. Q becomes: 



Energy conservation requires that (see Lahav 1991, Del Popolo & Gambera 1999): 



(8) 



(E) = (K) + (U G ) + (U A ) + <J7 L > = U GM + Uam + U hM 



(9) 



where the subscript "ta" stands for turn-around. The previous equation by means of Eq. (g) reads: 
iy^([/ G > + (2 + u){U A ) - u{U h ) = U OM + U AM + Uv M 

Eq. ^) and Eq. (|| or Eq. (|lo|), allows us to solve for K and r e «. 
In fact Eq. (bh can be written as: 



(10) 
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u + 1 



U G 



v + l\ GM 



1 + 2 



r c fs 
GAP 



r 3 3GM 



(11) 



where r e // = iprta, rt a being the radius of the turn-around epoch. Using the expressions: 



M = 4np b xl/3 



(12) 



£ = rt„/n,n = 2 |§S 1 Slgst;: 

° 

Inserting this last equation in Eq. (|ll|), I obtain: 



(13) 



(K) = 



(u+1) (GMH ) 2 ^Ql /a 



1 + 23. 



0€ 



2 2 1 /3^ 
Substituting G = in Eq. @, I finally get 

(v + 1) (GMH ) 2/3 fll /s 



GzQ.§H$M% Jo 



r3 3^0^ 



(14) 



(K) = 



10 2 

2^7^^ 



+ 3§n ^M§i -3 dr 3 fio-Ho 



(15) 



2 2V3^ 

where = ttW = 1 — Ho. Recalling the definitions of £, Eq. , and since at tun-around we have M = 47rptaT' 3 a /3, it is 
possible to express the product ifi£ as: 



r c s r ta r e ff I pb \ , N _i 
= (1 + «t a ) 



(see also Eq. (|24|) and Lilje (1992)). 

In order to get a connection between the kinetic energy and temperature, I utilize the usual relation: 
3/3MkT 



(16) 



(K) = 



2(im p 



(17) 



(AC), where k is the Boltzmann constant, /i = 0.59 is the mean molecular weight, m p the proton mass and f3 = kT J* m , 
being <r v the mass-weighted mean velocity dispersion of dark matter particles, and (5 = /3[1 + f(l//3 — l)Qb/Qm], where / is 
the fraction of the barionic matter in the hot gas, and Qb is the density parameter of the baryonic matter. 
Combining Eq. (|l|), Eq. (|l|) and Eq. @ I finally obtain: 



kT = i.58(^ + i)^— nl ( 1G M 



(1 + z ta ) 



10 2 
2 3 7T3 



1 



L , 2 
—z-dr 



A 



3t / """' J ff 2 noMf (1 + Zta ) 7 r " 3n ff 2 (l + 2ta) 3 
Eq. (|l8|) can be also equivalently written by using Lilje notation, in terms of r v j r : 

M 



m 3 



keV 



(18) 



kT 



0.94(1/ + 1)^2^) (^2. 



1 



/3 V?"vir / \ Pb,; 

lSr^^p^^ta 



3 1 / 

a* ' 



\10 15 MqH- 1 
" L 2 dr 2 A 



(1 + Zta) 



^// 2 f7oP? a ^L(l + ^ a )i r 3 3tf, 
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keV 



(19) 



8 A.Del Popolo et al. 





Fig. (3) Fig. (4) 

Figure 3. M-T relation predicted by the modified cluster continous-formation model. The solid line is the prediction of the quoted 
model, for Ha = 0, f2o = 1, shifted downwards, similarly to AC, to fit the FRB observational data in the massive end. The dotted and 
dashed line represents the prediction of the same model for L ^ 0, Q\ = 0, Qo = 1 an d L 7^ 0, S^a = and Cq = 0.3, respectively. 
Figure 4. Same as Fig. 3, but now the dotted line represents the prediction of the modified cluster continous-formation model for L y£ 0, 
A^0, (fi A = 0.7 and tt = 0.3). 



Using M = 3fp b A vir rl iT and M = ^pt a r 3 a , remembering that p b (z) = p [fio(l + z) 3 + (1 - O )] = poE(z) 2 and 



e2 ( z ) = n °n(t) )3 ( Br y an & Norman 1998), we see that: 



O 



-, 1/3 

1/3 



AX (1 + 2) (20) 



where 



A vir (z) ~ (18tt + 82y - 39y 2 ) (21) 
(Bryan & Norman 1998), being y = Q(z) — 1. These are accurate to 1% in the range Q(z) = 0.1 — 1. The value of Sl(z) is: 

" w n (i + *) 3 + (i-fio) ^ 

The previous result reduces to the standard one of late-formation approximation (see Eq. 8 of V2000) for v = L — A = 0. 

The value of r e ff, or the ratio ip = r e ff/r ta is obtained, similarly to Lahav et al (1991) by using Eq. ( |lo| ) which leads to 
the cubic equation: 

1 - *+ m 3 [n u + 2 V ) - * (2 + ^) - ("f + jH - (23) 

where 
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A Ar t 3 a 



2fi A 



( 



) 



-i 



(1+Zta) 



-3 



(24) 



n = 



47rGp t a 3GM 



Pta.b 



The previous equation reduces to the standard result = 1/2, for A = L = v = 0, flo = 1. 



The derivation of the previous relation is fundamentally based on the approximation of cluster formation with the evolution 
of a spherical top-hat density perturbation (Peebles 1993) and on the additional assumption that each cluster observed at 
a redshift z has just reached the moment of virialization. This last assumption is currently known as the late-formation 
approximation, which is a good one in a critical tto = 1, because for this value of Q massive clusters develop rapidly at all 
redshifts and the moment of virialization is always close to that of observation. In other terms for fio = 1, the accretion rate 
remains sufficiently high, and this implies that the clusters we actually observe attained their observed masses recently. In the 
flo < 1 case cluster formation is "shutting down" and it is necessary to take account of the differences between the moment 
of virialization and that of observation. The problem becomes worse going through Qo << 1: in fact in the late-formation 
approximation M v i r rises steadily since pbA v i r declines indefinitely, while we expect that the cluster formation is going to 
stop. As shown by V2000, to obtain the proper normalization and time evolution of the M-T relation, one has to account: 

a) for the continuous accretion of mass of clusters; 

b) for the non-zero density at r v i r , requiring a change in the virial theorem by including a surface pressure term. 

The M-T relation derived by means of a model of continuous accretion, as that of V2000, differs from the late-formation 
model in both normalization and time-dependent behavior. A comparison of the normalization predicted by the late-formation 
model with that predicted by simulations of Evrard, Metzler & Navarro (1996) shows that when Qo — 1 this normalization 
is only 4% below the empirical value, but it lies 20% below it for fio = 0.2. In the case of V2000 model and for a power-law 
spectrum, a comparison with the same simulations show that the temperature normalization of the n — —2 case deviates 
by less than 10% over the range 0.2 < Qq < 1 and by ~ 18% in the case n = — 1 (V2000). The normalization obtained by 
the V2000 model, even if it is more accurate than that given by the late formation, or that by AC which is in agreement 
with hydro-simulations, show a noteworthy discrepancy when compared with X-ray mass estimates (about 50% for the AC 
model; see also V2000). One possible source for differences in theoretical and observational normalizations may be due to the 
fact that P is different in the two cases because of systematic selection effects. For example, as shown by Bryan & Norman 
(1998), increasing the resolution of simulations there is an increase in the value of j3. So summarizing, for what concerns 
normalization, the continuous formation model gives more precise results than the late formation one, but in any case if we 
want to fit observations we need to shift the normalization (see AC). In any case, we have to stress that the continuous 
formation model must be considered as a noteworthy improvement on the late formation and is physically more consistent 
© 2000 RAS, MNRAS 000, §-?? 
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than the THM: as stressed by V2000, the top-hat model does not account properly for energy and mass accumulation during 
the early stages of cluster formation or the confining effects of matter that continues to fall in, both of which increase the 
temperature associated with a given mass (note that this last problem has been removed in the THM model). 

For the purpose of studying cluster evolution, the change in the M-T relation with z is more important than its normal- 
ization: although the differences in evolution of the M-T relation are of the order of 10% (V98), when it is introduced in the 
Press-Schechter model, they are amplified by the exponential term in the Press-Schechter formula. 

In the next section, I'll extend the V2000 model to take account of tidal interaction between clusters and the evolution 
of the M-T relation, comparing the results with the those of the top-hat model. 

3 REVISITING THE CONTINUOUS FORMATION MODEL 

In the previous section, I obtained the M-T relation for a universe with A ^ and L 7^ in the late-formation approximation. 
Some shortcomings of this approach has been summarized in the previous section (see also Viana & Liddle 1996; Kitayama 
& Suto 1996; Eke et al. 1996). The late-formation approximation is a good one for many purposes, but a better one can be 
obtained in the low-O limit. As can be found in the literature, there are two ways of improving the quoted model. One is 
to define a formation redshift z{ at which a cluster virializes and after the properties of observed clusters at z are obtained 
by integrating over the appropriate distribution of formation redshifts (Kitayama & Suto 1996; Viana & Liddle 1996). The 
second possibility is that described by V98, V2000. In this approach, the top-hat cluster formation model is substituted 
by a model of cluster formation from spherically symmetric perturbations with negative radial density gradients. The fact 
that clusters form gradually, and not instantaneously, is taken into account in the merging-halo formalism of Lacey & Cole 
(1993). In hierarchical models for structure formation, the growth of the largest clusters is quasi-continuous since these large 
objects are so rare that they almost never merge with another cluster of similar size (Lacey & Cole 1993). So, Lacey & Cole 
(1993) approach extends the Press-Schechter formalism by considering how clusters grow via accretion of smaller virialized 
objects. As shown by V2000, the mass grows like M oc w - 3 /(™+3) (Lacey & Cole 1993; V98; V2000). The virial energy of the 
cluster, —E, can be calculated by integrating the specific energy e, E — — J edM. The cluster temperature can be obtained 
remembering that it is proportional to Mg. Integrating Eq.(U), I get: 




(25) 



A particular shell will collapse if: 
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Fig. (5a) Fig. (5b) 

Figure 5a. Temperature evolution predicted by the top-hat M-T relation. The dotted line represents the prediction for f!o = 1 an d 
L = which coincides with the "classical" prediction, Tx oc (1 + z). The short-dashed line represents the case f!o = 0.2 and L = 0, 
coinciding with the late formation model described in V2000 (Eq. 8). The long-dashed line represents Qo = 1 and L ^ 0. Finally, the 
short-long-dashed line represents the case L ^ 0, A ^ (f2(> = 0.2, Ca = 0.8) and the solid line the case L ^ 0, A = 0, Co = 0.2. 
Figure 5b. Same as Fig. 5a but now the short-dashed line represents the case Qo = 0.3 and L = 0, while the short-long-dashed line 
represents the case L ^ 0, A / (Oo = 0.3, Qa = 0.7) and the solid line the case L ^ 0, A = 0, Qo = 0.3. 



GM 
To 



L a j. A 2 n 
M¥* + 6 ro= ° 



(26) 



(see V2000), and the shell reaches its maximum radius (turn-around radius, r ta ) at a time: 

dr 



; + GM rL^ dr+ A r2 
r Jr\ m*r a 6 



(27) 



Eq. ( |27[ ) can be written in an equivalent form (see Del Popolo & Gambera 1998, 1999; Bartlett & Silk 1993), as: 

dr 

tta 



y 2 dy 



(28) 



1 + 



<)?■> 



ATT 2 p 2 r}»(lH 2 Jl 



s/ 2 + (t£KV-[i+(^)£ 3 ] v 2 



The shell collapse at t c = 2tt&- In the case L = 0, Eq. ( p7[ ) can be analytically integrated and inverted to get e(t) 
1 (2ttGM\ 2/s _ 1 (2ttGM\ 2/3 (t n \ 2/3 1 (2nGM\ 2/3 



e(t) 



2 V t ) 



2 \ t a 



T (?) 



-J 



where we have defined in = 



2 V fa 

and x = 1 + (ii) 2 / 3 which is connected to mass by M — M x~ Zm/& (V2000) 



(29) 



H o (i-n -n A )! 

Integrating with respect to mass e(t), and diving by mass, one gets, as shown by V2000: 



E 
M 



f edM 



3m 



M 



2 

\ 3 2 

J Ms 


' 1 

h 


(?)*; 




m 





(30) 
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(6a) 

Figure 6a. Temperature evolution predicted by the modified cluster continuous-formation model in the case Ho = 0.2. The dotted line 
represents the "classical" prediction, Tx oc (1 + z). The short-dashed line represents the late-formation approximation as expressed by 

Eq. (8) in V2000, namely T x oc aV 3 [j^fcy] 1/3 (1 + z). The long-dashed line and the dot-dashed one plot Eq. @ for L = and 
n = —2, n = —1, respectively. The solid line plot the same equation for L ^ and n = — 1. 
Figure 6b. Same as Fig. 6a but for Qo = 0.3. 

In the case L 0, it is numerically possible to integrate and invert Eq. (^?]) finding the specific energy, which we can indicate 
with and finally: 

E 3 m 



-e h (t) 



(31) 



M 5 m - 1 ' 

An approximate relation for e^(t), in an Einstein-de Sitter Universe can be obtained as follows. It is possible to obtain the 
turn-around radius, rta by solving Eq. ( |2g| ) (with A = 0) for a given mass and a given epoch of interest. This is related to 
the binding energy of the shell enclosing the mass M by Eq. (^) with r = 0. In turn, the binding energy is uniquely given by 
the linear overdensity Si at some arbitrary early time. We may now use the linear theory to "propagate" the overdensity to 
the "chosen" time to find the linear overdensity at turn-around, S c . Using the relation between v and Si for the growing mode 
(Peebles 1980) in Eq. (^), with A = 0, at an early time, it is possible to connect e and 5 C (see Bartlett & Silk 1993)): 



e = -420^0*1 (oi/Oi) = - ( I 



(32) 



Using the formula for S c given in Del Popolo & Gambera (1999, Eq. 14), Del Popolo & Gambera (2000, Eq. 6), is possible to 
write: 



1 f2TvGM\ 2/3 r 1 f2TrGM\ 2/s 



2 V t ) 
which reduces to Eq. (tS2j) when L — *• 0. 



1 + 



GM 3 



Z/dr 



(33) 
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Eq. ( pS3| ) can also be written as: 

2^/3^/3^2/3 1 f r L2dr - 



1 /27rGM\ 2/3 , 

Defining M = M a;" 3m/5 and 
2 r WZ/& a [ r L 2 dr 



3 2 / 3 # 2 fiM 8 / 3 x - 1 



(34) 



F 



Eq. fl34h becomes: 



(35) 



/MX 



1 



1 



Integrating with respect mass, and dividing again by M, I get: 



(36) 



E _ fedM_ 3m (2nGy § 



M 



M 10(m - 1) V to / 



rn \ t ) M8/3 



(37) 



K(m,x) = (m- l)FxLerchPhi(x, l,3m/5 + l) - 



(m - 1) FLerchPhi(x, 1, 3m/5) 



where the LerchPhi function is defined as follows: 



(38) 



LerchPhi(z,a,v) = y — 
£ — ' (v 



(v + n) a 



(39) 



f\lf K = 0, Eq. (EtI) reduces to Eq. (10) of V2000. As stressed by V2000, some factors give rise to an higher value of E/M 
with respect the case of the late- formation value. The m/(m — I) value which accounts for the effect of early infall. The 1/m 
value in the square bracket of Eq. (^) which accounts for the cessation of cluster formation when t >> tn. Finally in Eq. 
( |37| ) a new term is present, which comes from the tidal interaction. In order to obtain an expression for the kinetic energy, 
starting from E/M, I use the virial theorem with the surface pressure term correction as in V2000 and Eq. (^). In this way, 
I finally get: 



. „ 2 um„ m /2ttG\3 2 
5 2/3 m - 1 V to / 



1 /to \ 3 K{m,x) 



M. 



8/3 



(40) 



where a = 2p ^ r p is the ratio between kinetic and total energy (V2000). Using the relation A v ir = j^r (see V2000), and 
in the early-time limit: (t << to), Eq. (fib]), reduces to: 

t This definition is valid for \z\ < 1. By analytic continuation, it is extended to the whole complex z-plane, for each value of a. If 
the coefficients of the series representation of a hypergeometric function are rational functions of the summation indices, then the 
hypergeometric function can be expressed as a linear sum of Lerch Phi functions. Reference: A. Erdelyi, 1953, Higher Transccndctal 
Functions, Volume 1, chapter 1, section 11. 
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2 m /imp 



* (f p b A vir ) 



1/3 



fcT = - 



5 m - 1 2/3 



(41) 



which, in the case n ~ —2, a ~ 2 is identical to the late-formation formula, described in V2000 (see their Eq. (8)). Normalizing 



where Ko(m, x) indicates that K(m,x) must be calculated assuming t = to 

Eq. (^) when compared to the result of V2000 (Eq. 17) shows an additional term, mass dependent. This means that, 
as in the case of the top-hat model, the M-T relation is no longer self-similar showing a break at the low mass end (see next 
section) . 

4 RESULTS 

The results of the models described in the previous sections are plotted in Fig. 1-6. 

Fig. 1 compares the M-T relation obtained by means of the THM model with the FRB data. The solid line is the 
prediction of Eq. (t^j), for J7a = 0, Qo = 1, shifted downwards, similarly to AC, to fit the observational data in the massive 
end. The reason of the shift is due to the fact that, as reported in the introduction, X-ray mass estimates lead to normalizations 
about 50% higher than hydro-simulations, (while mass estimates from galaxy velocity dispersion is consistent with simulation 
results (HMS)). The coincidence of the slope of the best fit to data and the theoretical model indicates that the process 
that the theoretical models does not take properly into account is probably happening at small scales (AC). The dotted line 
represents the prediction of Eq. (^) for i / 0, = 0, fio = 1. The plot shows a non-self-similar behavior. A break at 
Tx ~ 3keV is evident and in this low mass range the power index a ~ 1.8. It is interesting to note that the mass estimates 
using the observational /3-model lead to a power index 1.7 < a < 1.8. Originally this behavior was interpreted as an artifact 
of the /3-model (HMS), but the observation of a similar behavior for masses estimated from the resolved temperature profile 
(FRB) and confirmed by other studies (e.g., Xu, Jin & Wu (2001)), who found the break at Tx = 3 — 4keV. The quoted 
bent is interpreted in different ways. FRB suggests that it is due to different formation redshifts, but this conclusion is based 
upon the assumption that the temperature is constant after the formation time. This assumption is probably not true since 
clusters continue to accrete mass in time, even after cluster has formed. So as argued by Mathiesen (2001), using simulations 
results, the formation time might not be so important. Another interpretation is based on the assumption that the cluster 
medium is preheated in the early phase of cluster formation (see Xu, Jin & Wu 2001). In the AC study, the quoted break 



Eq. @ similarly to V2000, I get: 




(42) 
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was reproduced (justified), based on the scatter in the density field for a low density ACDM cosmology. The AC model is 
not able to distinguish between the effect of formation redshift from scatter in the initial energy of the cluster or its initial 
non-sphericity. In any case, non-sphericity introduces an asymmetric scatter dependent on mass for the M-T relation altering 
its slope at Tx ~ 3keV. As stressed by AC, the calculation of the effects of the non-spherical shape of the initial proto-cluster 
are not very rigorous and should be considered as an estimate of the actual corrections. 

In the model of this paper, the bent is entirely justified in terms of cluster tidal interaction with the neighboring ones, 
or in other terms it is strictly connected to the asphericity of clusters (see Del Popolo & Gambera 1999 for a discussion 
on the relation between angular momentum acquisition, asphericity and structure formation). Non-sphericity introduces an 
asymmetric bent, dependent on mass, in the M-T relation that gives rise to a different slope at the low mass end (T ~ 3keV): 
the lower the mass the larger the bent. A possible explanation for the reason of the origin of the bent is the following. The 
M-T relations found in this paper differ from previous models (e.g. V2000, AC), for the presence of a term depending on mass 
and angular momentum, L, originating from the gravitational interaction of the quadrupole moment of the protocluster with 
the tidal field of the matter of the neighboring protostructures. This last term changes the dependence of the temperature on 
the mass, M, in the M-T relation. The origin of the quoted term is due to: 

a) the change in the energetics of the collapse model and therefore of the turnaround epoch and the value of the threshold 
parameter S c that now is a function of the mass,M, produced by the introduction of another potential energy term ( A'j ) in 
iJq.Q. Angular momentum acquisition influence the value of S c which increases for peaks of low mass and remains unchanged 
for high mass peaks. Similarly to what is shown in Del Popolo & Gambera (1998, 1999), the threshold, <5 C , is a decreasing 
function of the mass, M. Since temperature, T, is directly proportional to the specific energy, e, of infalling matter and this in 
turn is directly proportional to the threshold parameter, S c , (see V2000), we have as a result that less massive clusters must 
be hotter while clusters in the high mass end have their temperature unchanged. 

b) the modification of the M-T relation produced by the alteration of the partition of energy in virial equilibrium. Moreover 
the partition of energy in virial equilibrium is altered and consequently also the M-T relation is modified. 

Coming back to Fig. 1, the dashed line represents the prediction of Eq. ( |l9| ) for Qa = and Qo = 0.3. Fig. 1 shows 

that in an Einstein-de Sitter universe the bent is under-estimated, while in a OCDM model with Qo = 0.3 it is overestimated. 
The comparison shows that the FRB data is able to rule out very low Qo models (< 0.3). 

Fig. 2 represents the same calculation of the previous figure but now Qa 0. The solid line is equivalent to the solid one 
in Fig. 1. The dashed line represents the case i / 0, A / 0, (Qo = 0.3, Qa = 0.7). It is evident that the effect of a non-zero 
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cosmological constant is that of reducing the effect of L, and that the ACDM cosmology, with Qo = 0.3, Qa = 0.7, is in better 
agreement with the observed bent than the Qo = 0.3 OCDM model. 

The M-T relation predicted by the modified continous formation model is plotted in Fig. 3 and Fig. 4. In Fig. 3, the solid 
line is the prediction of the quoted model, for Qa = 0, fio = 1, shifted downwards, similarly to AC and to the top-hat model 
(Fig. 1-2), to fit the FRB observational data in the massive end. The dotted and dashed line represents the prediction of the 
same model for L / 0, Qa = 0, Qo = 1 and L / 0, Qa = and Qo = 0.3, respectively. In agreement with the predictions 
of the top-hat model, the M-T relation shows a break at the low mass end (T ~ 3keV). The break is less "evident" in the 
case of a flat universe while much more evident in the case of an open one. A comparison with the FRB data shows, similarly 
to the top-hat model prediction, that in an Einstein-de Sitter universe the bent is under-estimated, while in a OCDM model 
with Qo = 0.3 it is overestimated. It is to be noted that the modified continuous formation model predicts a slightly larger 
bent than the top-hat model. 

Fig. 4 plots the M-T relation for the modified continuous formation model in a ACDM model. In this case, eL and then 
the M-T relation is calculated numerically by integrating and inverting Eq. (^). The solid line is the prediction for Qa = 0, 
Qo = 1, shifted downwards, as previously described and motivated, while the dotted line represents the prediction of the 
same model for L / 0, Qa = 0.7 and Qo = 0.3. The plot shows, similarly to the top-hat model, that the effect of a non-zero 
cosmological constant is that of reducing the effect of L, and that the ACDM cosmology (with Qa = 0.7 and Qo = 0.3) is in 
better agreement with the observed bent than the Qo = 0.3 OCDM model. The ACDM model predicts a slightly larger bent 
in the M-T relation when compared with the top-hat model. 

A comparison between the Fig. 1-2 and Fig. 3-4 shows that the prediction of the THM model and the VM model differs 
only for 1 — 2%. This little difference is possible because we have chosen the same normalization for the two models. 

The study of the evolution of the M-T relation predicted by the modified top-hat model, for a fixed value of M v i r , is 
plotted in Fig. 5a-b. 

In Fig. 5a, the dotted line represents the prediction for Q — 1 and L — which coincides with the "classical" prediction, 
Tx oc (1 + z). The short-dashed line represents the case Qo = 0.2 and L = 0, which coincides with Tx oc AV r 3 
given in V2000 (Eq. 8). As known, (see V98), decreasing Q slows down the evolution. The long-dashed line represents Qo = 1 
and Comparing this result to the first one Qo — 1 and L — shows that the effect of L in slowing down the evolution of 

Tx is larger than reducing Qo of a noteworthy value (0.8, see short-dashed line). Finally, the short- long-dashed line represents 
the case L 7^ 0, A 7^ 0, Q — 0.2, Qa = 0.8 and the solid line the case L 7^ 0, A = Qo = 0.2. As expected, reducing the value 
of Q and taking account of tidal interaction of clusters (L 7^ 0) produces a larger effect in the slowing down of evolution, while 
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the effect of a non-zero cosmological constant has an opposite effect to that of L. Fig. 5b plots the same calculations of Fig. 
5a but for the case £7q = 0.3. As expected, an increase in produce a more rapid evolution of Tx- 

Fig. 6a, is based upon the second model (modified continuous formation) and calculates the evolution for Qo = 0.2. 
The dotted line represents the "classical" prediction, Tx oc (1 + z). This relation predicts the largest evolution for a given 
Mvir. The short-dashed line represents the late-formation approximation as expressed by Eq. (8) in V2000, namely Tx oc 
(1 + z). In this case, temperature evolution is more modest than the previous one. The long-dashed line and 
the dot-dashed one, plot Eq. (^) for L = and n = —2, n = — 1, respectively. This last two lines show that, (see also V98), 
the continous formation model produces an even more modest evolution of temperature with respect to the late-formation 
approximation. The solid line plot Eq. ( |4o| ) for L / and n = — 1. The effect of the angular momentum is that of furtherly 
increase the slowing down effect of the temperature evolution. Fig. 6b, plots the same calculations of Fig. 6a for fio = 0.3. 

A comparison between the Fig. 5 and Fig. 6 shows that the prediction, for what concerns the time evolution, of the THM 
model and the VM model differs < 10% (for L / 0, A = 0, fi = 0.2). 



5 CONCLUSIONS 

In this paper, I used two different models in order to study the M-T relation and its time evolution. The first one is based 
on a modification of the top-hat model in order to take account of an external pressure term in the virial theorem and of 
angular momentum acquisition by protostructures (see Del Popolo & Gambera 1999) , the second one is a modification of the 
continuous formation model of V2000, again to take account of angular momentum acquisition by protostructures. 
The main results of the paper can be summarized as follows: 

1) the effect of angular momentum acquisition by protostructures is two fold: 

a) The M-T relation is no longer self-similar: a break in the low mass end (T ~ 3 — 4kcV) of the M-T relation is present. 
The behavior of the M-T relation is as usual, M oc T 3 ' 8 , at the high mass end, and M oc T 7 , with a value of 7 > 3/2 in 
dependence of the chosen cosmology. Larger values of 7 are related to open cosmologies, while ACDM cosmologies give results 
of the slope intermediate between the flat case and the open case. The FRB data are best fitted by ACDM models and by 
CDM models with Q > 0.3. 

2) The evolution of the M-T relation is more rapid in models with L = 0. The effect of a non-zero cosmological constant is 
that of slightly increase the evolution of the M-T relation with respect open models with the same value of £Iq This is due 
to the fact that in a flat Q,a > universe, the universe's density remains close to the critical value later in time, promoting 
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perturbation growth at lower redshift. The evolution is more rapid for larger values (in absolute value) of the spectral index, 
ri. 

3) The top-hat model gives comparable results to the continuous formation model, in the range of z and Qo considered. For 
sake of precision, the top-hat model has a lower value of normalization. A comparison of the normalization predicted by the 
late-formation model with that predicted by simulations of (Evrard, Metzler & Navarro 1996) shows that when Qq = 1 this 
normalization is only 4% below the empirical value, but it lies 20% below it for Qq = 0.2. In the case of V2000 model and for 
a power-law spectrum, a comparison with the same simulations show that the temperature normalization of the n = —2 case 
deviates by less than 10% over the range 0.2 < f2o < 1 and by ~ 18% in the case n = — 1 (V2000). Given the same value of 
normalization for both models, a comparison between the Fig. 1-2 and Fig. 3-4 shows that the prediction of the THM model 
and the VM model differs only for 1 — 2%. A comparison between the Fig. 5 and Fig. 6 shows that the prediction, for what 
concerns the time evolution, of the THM model and the VM model differs < 10% (for L 7^ 0, A = 0, Qo = 0.2). I want to stress 
that even if in the parameter space considered, the THM model gives similar results to those of the continuous formation 
model, one cannot forget that it does not account properly for energy and mass accumulation during the early stages of cluster 
formation. In any case the central aim of this paper is to study the M-T relation and its evolution, the comparison between 
the THM and VM is only a secondary point. 
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